Crossing the Phantom Divide: Theoretical Imphcations and Observational Status 
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If the dark energy equation of state parameter w{z) crosses the phantom divide hne w — —1 (or 

equivalently if the expression '^^^^^'^^^ — 3r2omJ^o(l + ^)^ changes sign) at recent redshifts, then there 
are two possible cosmological imphcations: Either the dark energy consists of multiple components 
with at least one non-canonical phantom component or general relativity needs to be extended 
to a more general theory on cosmological scales. The former possibility requires the existence 
of a phantom component which has been shown to suffer from serious theoretical problems and 
instabilities. Therefore, the latter possibility is the simplest realistic theoretical framework in which 
such a crossing can be realized. After providing a pedagogical description of various dark energy 
observational probes, we use a set of such probes (including the Gold Snia sample, the first year 
SNLS dataset, the 3-year WMAP CMB shift parameter, the SDSS baryon acoustic oscillations 
peak (BAO), the X-ray gas mass fraction in clusters and the linear growth rate of perturbations 
at 2 = 0.15 as obtained from the 2dF galaxy redshift survey) to investigate the priors required for 
cosmological observations to favor crossing of the phantom divide. We find that a low fiom prior 
(0.2 < f2om < 0.25) leads, for most observational probes (except of the SNLS data), to an increased 
probability (mild trend) for phantom divide crossing. An interesting degeneracy of the ISW effect 
in the CMB perturbation spectrum is also pointed out. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

The assumption of large scale homogeneity and 
isotropy of the universe combined with the assumption 
that general relativity is the correct theory on cosmolog- 
ical scales leads to the Friedman equation which in a flat 
universe takes the form 



[l6j for recent reviews) needed to induce the accelerat- 
ing expansion are well described by its equation of state 
w{z) = ^^1^1 which enters in the second Friedman equa- 
tion as 
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where a{t) is the scale factor of the universe and p 
its average energy density. Both sides of this equation 
can be observationally probed directly: The left side 
using mainly geometrical methods (measuring the lu- 
minosity and angular diameter distances dt(2;)[l|, 0, Q 
and dA{zM 1 1 

with standard candles and stan- 
dard rulers) showing an accelerating expansion at re- 
cent redshifts and the matter - radiation density part 
of the right side using dynamical and other methods 
(cosmic microwave background jTi], large scale structure 
observations 0, lensingj^ etc). These observations have 
indicated [13] that the two sides of the Friedman equa- 
tion (jl.ip can not be equal if p{a) — Pm{a) ~ 
even if a non-zero curvature is assumed. There are 
two possible resolutions to this puzzle: Either modify 
the right side of the Friedman equation Ijl.ip introduc- 
ing a new form of 'dark' energy ideal fluid component 
{p{a) = Pm{a) + Pxio)) with suitable evolution in order 
to restore the equality or modify both sides by changing 
the way energy density affects geometry thus modifying 
the Einstein equations. 

In the first class of approaches the required gravita- 
tional properties of dark energy (see [ll|, [13, fl3l Il4 [Tsl . 



implying that a negative pressure {w < —1/3) is nec- 
essary in order to induce accelerating expansion. The 
simplest viable example of dark energy is the cosmologi- 
cal constant [l^, [13, [l3l {w = — 1 ) . This example however 
even though consistent with present data lacks physical 
motivation. Questions like 'What is the origin of the 
cosmological constant?' or 'Why is the cosmological con- 
stant 10^^° times smaller than its natural scale so that 
it starts dominating at recent cosmolo gica l times (coin- 
cidence problem)?' remain unanswered [l9j. Attempts to 
replace the cosmological constant by a dynamical scalar 
field (quintessence [20,, .21, 22]) have created a new prob- 
lem regarding the initial conditions of quintessence which 
even though can be resolved in particular cases (tracker 
quintessence), can not answer the above questions in a 
satisfactory way. 

The parameter w{z) determines not only the gravita- 
tional properties of dark energy but also its evolution. 
This evolution is easily obtained from the energy mo- 
mentum conservation 



d{pxa^) = -pxd{a^) 



(1.3) 



which leads to 
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Therefore the determination of w{z) is equivalent to that 
of Px{z) which in turn is equivalent to the observed H[z) 
from the Friedman equation (jl.ip expressed as 



H{z) = Ho[nonA^ + zf + noxe 



3/0 ii^li+^'l^'))! 



(1.5) 



Thus, knowledge of ilom and H{z) sufRces to determine 
u>(z) which is obtained from equation ()1.5p as [13] 
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In the second class of approaches the Einstein equa- 
tions get modified and the new equations combined with 
the assumption of homogeneity and isotropy lead to a 
generalized Friedman equation of the form 



(1.7) 



where / and g are appropriate functions determined by 
the modified gravity theory [H, [H, HI, 113, [M IS! • In this 
class of models, the parameter w{z) can also be defined 
from equation (jl.6p but it can not be interpreted as ^ 
of a perfect fluid. 

The simplest (but quite general) examples of modi- 
fied gravity theories are scalar tensor theories [23. IsoL IbiI. 
[3^ . [33} where the Newton's constant G is promoted to a 
function of a field SttG — > whose dynamics at 

the Lagrangian level is determined by a potential U{^). 
Assuming homogeneity and isotropy, the modified Fried- 
man equation in these theories take the form 



+ U -3HF) (1.8) 



This equation reduces to a regular minimally coupled 
scalar field dark energy (quintessence) in the general rel- 
ativity limit of a constant F = . The positive nature 
of the kinetic term $^ however implies that certain types 
of behaviors of H may not be reproducible without in- 
voking a time-dependent F. These types of behavior of 
H{z) which include a w{z) crossing the Phantom Divide 
Line (PDL) w = —1 are potential signatures of extended 
gravity theories and will be discussed in the next section. 

To identify this type of signatures, a detailed form of 
the observed H(z) is required which may be obtained 
by a combination of multiple dark energy probes. Such 
probes may be divided in two classes [IJl according to the 
methods used to obtain H{z). 

• Geometric methods probe the large scale geometry 
of space-time directly through the redshift depen- 
dence of cosmological distances {dL{z) or dyi(z)). 
They thus determine H{z) independent of the va- 
lidity of Einstein equations. 

• Dynamical methods determine H{z) by measuring 
the evolution of energy density (background or per- 
turbations) and using a gravity theory to relate 



them with geometry ie with H{z). These methods 
rely on knowledge of the dynamical equations that 
connect geometry with energy and may therefore 
be used in combination with geometric methods to 
test these dynamical equations. 

Examples of geometric probes include 

1. The measured supernova distance redshift relation 
t^L(-z)[il, 0,131 which for a fiat universe, is connected 
to H{z) as 



dL{z) ^(l + z) 
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2. The measured^, HI] angular diameter distance 
dA{zrec) to the sound horizon rs{zrec) at recom- 
bination 



dAi^rec) — 
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3. The scale of the sound horizon measured at more 
recent redshifts {zbao) through large scale struc- 
ture redshift survey correlation functions Q 
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1/3 



4. The cluster gas mass fraction defined as 
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assumed to be constant for all clusters and pro- 
portional to probes the angular diameter 
distance [i, H, [HI to each cluster as discussed in 
detail in section III. 

An example of a dynamical probe of geometry is the mea- 
sured linear growth factor of the matter density pertur- 
bations D{a) defined as 
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The measurements of F){cl) can be made by several meth- 
ods including the redshift distortion factor in redshift 
surveys ^3 J weak lensing[38[, number counts of galaxy 
clustersjUl, Integrated Sachs- Wolfe (IS W)effect fw] and 
large scale structure power spectrum[4l|, [l^. The theo- 
retical prediction of the evolution of -D(a) on sub-Hubble 
scales is obtained from the Euler and matter stress energy 
conservation equations asfH, HI] 
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FIG. 1: The growth rate in the (flat and Qom = 0.26) ACDM 
(continuous line), Scalar- Tensor (long-dashed line) and DGP 
(short-dashed line) cases and a datapoint from the 2dFGRS 
g{a = 0.15) = 0.51 ± 0.11. The ACDM model is consistent 
with the current growth rate observations. 



with initial conditions D{a) ~ a for a ~ 0. In equa- 
tion p.l4p we have ignored anisotropic stresses and dark 
energy perturbations |l^ which are expected to have a 
small effect on sub-Hubble scales. The last term of equa- 
tion (frT4|) emerges by connecting the metric perturba- 
tion with the matter density perturbations. It therefore 
depends on the particular form of the dynamical equa- 
tions of the gravity theory considered. This dependence 
is expressed through the function f{k, a) which in the 
case of general relativity is unity (f{k,a) = 1) while in 
extended gravity theories it can take values different from 
one which can even depend on the scale k [i^ . 
For example for scalar-tensor theories we have 



f{k,a) 



Geff{a) 




(1.15) 



where Geff{a) is the effective Newton's constant when 
the scale factor is a and a — 1 corresponds to the present 
value of the scale factor while m is the mass of the scalar 
field $ inducing a Yukawa cutoff to the gravitational field. 
In what follows (Fig. 1) we will use a simple ansatz for 
the effective Newton's constant 
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and neglect the mass m of the scalar field. Alternatively, 
for the DGP modelfHEl we havelH 



with 
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where Vc is the crossover scale beyond which the gravita- 
tional force follows the 5-dimensional l/r^ behavior 



(1.19) 



and H{a) for the DGP model is defined through [41 
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Solving for Hdgp{o), assuming flatness and in the case 
when we only have matter on the brane we get 



H 



DGP 



(a) 



where now 



(1.21) 



fl,^ = -(l-f7o„0' (1.22) 

The detection of an f{k, a) ^ 1 from equation (|1.14[) 
would therefore be a 'smoking gun' signature of extended 
gravity theories. Such a detection could be made for ex- 
ample by using the form of H(z) obtained from geometric 
tests in equation (|1.14p . solving for D{a) in the context 
of general relativity (/(fc, a) — 1) and comparing with 
the observed D{a) at various redshifts. If a statistically 
significant difference is found between the observed D{a) 
and one predicted in the context of general relativity then 
this could be interpreted as evidence for extensions of 
general relativity. 

There is currently an observational estimate of the 
growth rate defined as 



5(a) 



aD'{a) 
D{a) 



(1.23) 



at a redshift z = i-1 = 0.15 from the 2dFGRS [1, [13, H 
It is 



g(z = 0.15) = 0.51 ±0.11 



(1.24) 



and its derivation used the redshift distortion factor (see 
section III E) . There will soon be better estimates coming 
from the SDSS^41j|. In Fig. 1 we show the growth rate for 
the following flat models: a best fit (fiom = 0.26) ACDM 
model, a DGP (l^om = 0.22) model and a Scalar- Tensor 
model (^ = —0.2 in equation p.l6p and fiom = 0.26). As 
demonstrated in Fig. 1 however, the large errorbars in 
the currently available datapoint of equation (|1.24p can 
not distinguish the best fit ACDM model (flat f^A = 0.74) 
from other competing models. We therefore conclude 
that several different models including modified gravity 
models as well as ACDM are in agreement with current 
data. 

Given the current uncertainties in the growth rate ob- 
servations, it is important to uncover potential signatures 
of extended gravity theories in geometric observational 
methods. As discussed in the next section, a minimally 
coupled scalar field dark energy can not reproduce a w{z) 
crossing the PDL (w = —1) for any scalar field potential 
(see however Ref [47[ for an interesting case). The sim- 
plest realistic theoretical model which can reproduce such 
crossing is a scalar-tensor extension of general relativity. 
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Therefore, if such crossing is confirmed by geometric ob- 
servations then this could be interpreted as an indication 
for extended gravity theories. 

The structure of this paper is the following: In the next 
section we discuss the theoretical consequences emerging 
if the crossing of the PDL is confirmed by future observa- 
tions. In section III we fit the parameterization [isl. [49l| 



w{z) — Wq + Wi 



l + z 



(1.25) 



to current data from various dark energy observa- 
tional probes including the Gold Snia sample 0, the 
first year SNLS datasetg, the 3-year WMAP CMB 
shift parameter 0, [s^, the SDSS baryon acoustic peak 
(BAG) 3, the X-ray gas mass fraction in clusters [51, |6| 
and the linear growth rate of perturbations at z = 0.15 
as obtained from the 2dF galaxy redshift survey [1, [13], to 
investigate which of these probes currently favor crossing 
of the PDL and what are the priors required. Finally 
in section IV we conclude, summarize and outline future 
prospects of this work. 



II. THEORETICAL IMPLICATIONS OF PDL 
CROSSING 



The simplest class of physically motivated models gen- 
eralizing the cosmological constant and producing accel- 
erated universe expansion are based on a siniple evolving 
homogeneous scalar field $ (quintessence ^ [H HI]) 
minimally coupled to gravity whose dynamics is deter- 
mined by a potential ?7($) and has a canonical kinetic 
term It is easy to see however that this class of mod- 
els can not reproduce a 'w{z) crossing the PDL w = — 1 
for any potential U{^). Indeed, the equation of state w 
for such models takes the form 



w{z) 



px 

PX 



1$2 
2^ 



i$2 



(2.1) 



which approaches w = — 1 in the limit of a small kinetic 
term but does not cross the PDL as long as <I>^ > 
for any sign of J7($). This result has been general- 
ized by Vikmanfs^ who showed that any minimally cou- 
pled scalar field with a generalized kinetic term also can 
not cross the PDL through a stable trajectory (see Refs 
[Sll [5^ for an interesting exception involving however an 
arbitrary change of the kinetic term sign) . 

A simple 'no go' theoremfsl, [sj] can also be obtained 
for a general perfect barotropic fluid with conserved en- 
ergy momentum tensor whose equation of state is of the 
form p = p{p) . Indeed energy conservation implies that 



-3-(p + p(p)) 
a 



From this equation it follows that p ^ as p{p) 
as w = ^ ^ — 1. By differentiating equation 



(2.2) 



-p or 
with 



respect to time it is easy to see that p ^ Q as w — 1 
and similarly for all time derivatives of the density p ie 



lim 

w — * — 1 







(2.3) 



for all n. Therefore, any fluid with equation of state of 
the form p — p{p) can not cross the PDL. Instead it 
asymptotically approaches a constant energy density p 
thus mimicking a cosmological constant. It should be 
stressed however that many scalar field models do not 
behave like barotropic [p = p{p)) fluids due to the de- 
pendence of p and p on two independent variables $ and 
$ and therefore this proof is not equivalent to the proof 
that a minimally coupled scalar fleld can not cross the 
PDL. 

Since the minimal theoretical approaches are unable to 
reproduce a PDL crossing we must turn to 'non-minimal' 
models to achieve such crossing. There are two main 
classes of such 'non-minimal' approaches: Either consider 
multiple component dark energy [54,, i55|, i56i, i57|, i58i, iSl, 
113, [fill, m, mill Hi (notice that higher derivative mod- 
els and many non-barotropic fluids are effectively multi- 
component) with at least one phantom degree of freedom 
(eg scalar field with negative kinetic ener gy) or consider 
extensions of general relativity [23 . [25l . [2^ l27l. [28[ . [sol [H, 
[li,[63, 68, 69, 7|[lll[2i,[2l,[2l,[zl[2fl[2l- The phantom 
degrees of freedom of the former approach expressed by 
phantom fields are plagued by catastrophic UV instabil- 
ities since their energy is unbounded from below and al- 
lows vacuum decay through the production of high energy 
real particles and negative energy ghosts [78,, X9, 80, 81|]. 
Gn the other hand the general relativity extensions of 
the latter approach are severely constrained by local so- 
lar system and by cosmological observations but are well 
motivated theoretically. 

To illustrate the multiple dark energy component ap- 
proach to PDL crossing with scalar fields let's consider a 
set of two coupled real scalar fields (quintessence $1 -I- 
phantom $2 '■ quintom dark energy [57] ) with Lagrangian 

^ = \{d^.^if - \{d,^2f - C/($i, $ 2) (2.4) 

The effective pressure and energy density for a homoge- 
neous system is 



leading to the equation of state parameter 



w 



$2 



$2 



2C/($i,$2) 



$2 



$2 



2t/($i,$2) 



(2.5) 
(2.6) 

(2.7) 



<i>2 changes sign. 



which crosses the PDL line when $^ 
This can easily be achieved with appropriate potentials 
and initial conditions. 
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The same approach may be illustrated by conside ring: 
a mixture of multiple perfect fluids [H, [Sa, [lO, [H, [s^l 
instead of scalar fields. Consider a mixture of two non- 
interacting fluids and {p2,P2,W2) with sepa- 
rately conserved energies and constant equation of state 
parameters wi, W2 with wi > —1, W2 < —1. The equa- 
tion of state parameter w of the mixture is [s^ l 



w = — = 



Pi +P2 
Pi + P2 



wipoia 



3(l+tui) _|_ 



W2Pma 



Poia 



(2.8) 

which interpolates between w = wi {a << 1) and w = W2 
{a » 1) thus crossing the PDL u; = — 1. 

The requirement of phantom degrees of freedom which 
are plagued with several unattractive features [tsI. FtqI. [soI. 
[mIj (instabilities and lack of reahstic prototypes) 

makes this class of approaches theoretically unattractive 
compared to the second class which is based on extensions 
of general relativity. Even though the latter are severely 
constrained observationally they are strongly motivated 
from the theoretical viewpoint. First it is clear that gen- 
eral relativity is an incomplete theory because it does 
not contain quantum mechanics and also it is plagued by 
singularities. Second, all theories that attempt to quan- 
tize gravity and/or unify it with other interactions re- 
quire modifications of general relativity. For example, 
the string theory dilaton [ll|, [s^ field can be understood 
as the scalar field of an effective scalar-tensor theory in 4- 
dimensionsfSG'l. Similarly, all Kaluza-Klein type theories 
which attempt to unify gravity with other interactions 
by utilizing compactified extra dimensions are described 
as effective scalar-tensor theories at low energies [13, [11]. 

The simplest but very general (given its simplicity) ex- 
tension of general relativity is expressed through scalar- 
tensor theories. In these theories Newton's constant ob- 
tains dynamical properties expressed through the poten- 
tial F{^). The dynamics are determined by the La- 
grangian density [2J, [sO] 

F(^) 1 

jC = R-- g'^-'df.^d,^ - [/($) + Cni^Pm-, ffM.] 

. (2.9) 

where £m[V'm; duv] represents matter fields approximated 
by a pressureless perfect fiuid. The function F{^) is ob- 
servationally constrained as follows: 

• F($) > so that gravitons carry positive 
energvfsoj. 

• ^ < 10~^ from solar system observations [s^. 

Assuming a homogeneous $ and varying the action cor- 
responding to (|2.9p in a background of a flat FRW met- 
ric, we find the coupled system of generalized Friedman 
equations 



3FH^ = p+^^^ -3HF + U 



-2FH 



p + p + <^^+F-HF 



(2.10) 
(2.11) 



where we have assumed the presence of a perfect fluid 
(p = Pm,P — 0) playing the role of matter fields. Ex- 
pressing in terms of redshift and eliminating the potential 
U from equations ([TTU)) . ^(TH}, we find [laii] 



$'2 ^-F" - 



(InH)' 



Ho 



1 + z 



(InH)' 
F' '-F- 



1 + z 



-3(1 + z)r!om ( ) ^^0 > 



(2.12) 



where ' denotes derivative with respect to redshift and Fq 
is set to 1 in units of g^^^ and corresponds to the present 
value of F. In the limit of minimal coupling {F = 1) 
equation (j2.12ll becomes 



$'2 (H) 



;^^^-3(l-hz)r!om(-^)' >0 (2.13) 



Using equation (|1.6p it is easy to show that the inequality 
(|2.13p is equivalent to w{z) > —1 which confirms the 
expectation that a minimally coupled scalar field is not 
consistent with PDL crossing. This constraint however 
is relaxed due to the allowed dynamical evolution of F in 
equation (|2.12p which may be written as 



^'\f,h) = ^'„Uh) 



A(F, H)>0 



(2.14) 



where 



A{F, H) = F" 



{InH)' 



1 



F' 



(InH)', 
2 \ ' (l-F) 



Therefore, the constraint 



$'2 

nic\ 
'2 



l + z 

(2.15) 

„{H) > is replaced in 
scalar-tensor theories by ^^^{H) > A(F, H) and by 
choosing an F such that A(F, H) < a PDL crossing 
can be achieved. 

The question to be addressed is the following: 'Do ob- 
servational constraints on F allow A{F, H) < 0?' The 
answer to that question is positive as may be easily seen 
by considering the present time {z = 0) when the ob- 
servational constraints [qO] (and especially the solar sys- 
tem tests [1^) are most stringent requiring F'{z) ~ but 
placing no constraints on F"{z = 0). Setting z = 0, 
F = Fq ^ I, F' = Q in equation fTT^ we obtain 



w{z = 0) + 1 > 



0) > 3nom + F"{z = 0) 
F"{z = 0) 



3(1 - r!om) 



(2.16) 



where we have used equation (|1.6p . Since there is no 
observational constraint on F"{z ~ 0) we can clearly 
pick F"{z = 0) < allowing for w{z = 0) < — 1 and 
crossing the PDL. 

There is a simple physical interpretation of the behav- 
ior required by F in order to cross the PDL and lead 
to superacceleration {w < —1). Equation (|2.15p implies 
that A{H,F) < (and therefore w < —1) is favored for 



6 



F" < 0, F' < and F > 1. Since F ~ G~jj this behavior 
imphes an effective Newton's constant that increases with 
redshift (G"jy > 0, G'^ff > 0) and therefore decreases 

with time (Ge// < 0, Ge// < 0, Ge//(i) > Ge//(to)) 
thus 'helping' the accelerating expansion induced by the 
potential U. This type of behavior was verified in a spe- 
cific reconstruction example of F{z), U{z) in Ref. |33l | 
corresponding to an H{z) — w{z) that crosses the PDL. 
It should be stressed however that any ansatz of Ge// 
that is monotonic function of time, eg Ge//(i) ^ 
where a > 0, can lead to very tight constraints on the 
present values of G^j^ and G'^jf using nucleosynthesis 
constraints [91] (but see also [o^ for a way to relax such 
constraints). 

To summarize we have discussed two broad classes of 
models which can lead to crossing of the PDL: multi- 
component dark energy with phantom components and 
extensions of general relativity. This is not an ex- 
haustive classification but many of the left out models 
can be incorporated in the above classes. For exam- 
ple models of coupled quintessence where matter density 
[631 is explicitly coupled to the scalar field causing ac- 
celeration may be shown to be conformally equivalent 
to scalar-tensor theories. Similarly, braneworld mod- 
els [26, ,27,, ed, m. III, [9I, [Hi can be classified as ex- 
tended gravity theories. On the other hand in the multi- 
component class we could include models with a complex 
equation of statefU, [95| . fermion models, scalar fields 
with higher derivatives |64l| and vector field models [B^I- 
Of the two classes of theories discussed, the class of ex- 
tended gravity theories is strongly motivated theoreti- 
cally in contrast to the multicomponcnt class. There- 
fore, despite the observational constraints, extensions of 
general relativity is the prime candidate class of theories 
consistent with PDL crossing. 

It is thus important to address the following questions: 

• Do current cosmological data support a crossing of 
the PDL? 

• What are the optimal observational strategies to 
confirm or exclude a PDL crossing with future ob- 
servations? 

The main goal of the next section is to address the 
first question using a broad sample of cosmological data. 
Previous studies [Oq have indicated that the Snia data 
do not show a common trend regarding the crossing 
of the PDL. In particular, the analyses of the Gold 
SnIa dataset have indicated that a w{z) that crosses 
the PDL is preferred over the ACDM param eter val- 
ues at a level of almost 2a [13, [H, ^ [To^. On the 
other hand the SNLS da taset does not show such a trend 
^01, 102, 103, [io3, [lOil and the ACDM parameter val- 
ues are well within la from the best fit parameter values. 
This raises the question as to what is the trend favored by 
the other than SnIa cosmological data. Do they favor the 
trend of PDL crossing indicated by the Gold dataset or 
do they favor a constant w ~ — 1 as indicated by SNLS? 



The main goal of the next section will be to address this 
question using the following cosmological data: 

• The location of the CM B perturbatio n spectrum 
peaks0 (shift parameter [loi, [lol, [loi|). 

• The measurement of the CMB acoustic scale at 
zbao — 0.35 as indicated by the large scale struc- 
ture correlation function of SDSS redshift survey 
(Baryon Acoustic Oscillations (BAO) peak) [4]. 

• The gas mass fraction of galaxy clusters fgas — 

inferred from X-ray observations [1, Q , which 
depends on the angular diameter distance cIa to the 
cluster as fgas ^ f^A^ ■ 

• The linear growth factor D{a) (cq. ll.lSp as de- 
termined by the 2dF galaxy redshift survey [1, W?\ 
which can constrain H{z) by using eq. ()1.14p in the 
context of general relativity. 

Using the above cosmological data we find the best fit 
form of w{z) with la errors in the context of the GPL 
parameterization [isL l49l | 

w(a) = WQ + wi{l — a) (2-17) 

w{z) = wo + wi-^ (2.18) 
1 -|— z 

We then compare this form with the corresponding best 
fits of both the Gold and the SNLS datasets to see which 
of the two trends is favored. The comparison is made for 
priors of fiom in the range 0.2 — 0.3. 



III. OBSERVATIONAL OF PDL CROSSING 

A. The SnIa datasets 

The two most reliable and robust SnIa datasets ex- 
isting at present are the Gold dataset [1] and the Su- 
pernova Legacy Survey (SNLS) [1] dataset. The Gold 
dataset compiled by Riess et. al. is a set of supernova 
data from various sources analyzed in a consistent and 
robust manner with reduced calibration errors arising 
from systematics. It contains 143 points from previously 
published data plus 14 points with z > 1 discovered re- 
cently with the HST. The SNLS is a 5-year survey of 
SnIa with z < 1. It has recently @ released the first 
year dataset. The SNLS has adopted a more efficient 
SnIa search strategy involving a 'rolling search' mode 
where a given field is observed every third or fourth night 
using a single imaging instrument thus reducing pho- 
tometric systematic uncertainties. The published first 
year SNLS dataset consists of 44 previously published 
nearby SnIa with 0.015 < z < 0.125 plus 73 distant SnIa 
(0.15 < z < 1) discovered by SNLS two of which are out- 
liers and are not used in the analysis. The fact that in 
the two datasets a set of low-z SnIa is common to both 
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samples could only lead to minor common systematics 
due to low redshift. 

The above observations provide the apparent mag- 
nitude m{z) of the supernovae at peak brightness af- 
ter implementing correction for galactic extinction, K- 
correction and light curve width-luminosity correction. 
The resulting apparent magnitude m(z) is related to the 
luminosity distance Dl(z) through 

mth{z) = M{M,Ho) + 5logio{DL{z)) (3.1) 

where in a flat cosmological model 

DLiz) = {l + z) r dz' , ^° ^ (3.2) 

is the Hubble free luminosity distance {Hod^), ai, ...,a„ 
are theoretical model parameters and M is the magnitude 
zero point offset and depends on the absolute magnitude 
M and on the present Hubble parameter Hq as 

M ^ M + 5logw{-^) + 25 = 
Mpc 

= M ~ 5logioh + A2.38 (3.3) 

The parameter M is the absolute magnitude which is 
assumed to be constant after the above mentioned cor- 
rections have been implemented in m{z). 

The data points of the Gold dataset are given after 
the corrections have been implemented, in terms of the 
distance modulus 

fi'^,^{zi) = m^,^{z,) - M (3.4) 

The SNLS dataset however also presents for each point, 
the stretch factor s used to calibrate the absolute magni- 
tude and the rest frame color parameter c which mainly 
measures host galaxy extinction by dust. Thus, the dis- 
tance modulus in this case depends apart from the abso- 
lute magnitude M, on two additional parameters a and 
f3 defined from 

= ^Ts'-'iz,) - M + «(., - 1) - /3c, (3.5) 

which are fit along with the theoretical parameters using 
a recursive procedure. 

The theoretical model parameters are determined by 
minimizing the quantity 

XSnIa[^Om,Wo,Wi) ^ y ^ (3.6) 

where N = 115 for SNLS {^snls)^ ^ = 157 for the Gold 
dataset (jdooid)^ ^^li^ '^Lt and , are the errors due to 
flux uncertainties, intrinsic dispersion of Snia absolute 
magnitude and peculiar velocity dispersion respectively. 
These errors are assumed to be gaussian and uncorre- 
lated. The theoretical distance modulus is defined as 

^J■thiz^) = rnth{zi) - M = blogio{D l{z)) + /io (3.7) 




0.5 1 1.5 1.75 0.5 1 1.5 1.75 



FIG. 2: The best fit form of w{z) for the Gold and SNLS 
datasets for a prior of f2om = 0.24 along with the la errors 
(shaded region). 

where 

Aio = 42.38 - hlogwh (3.8) 

and ^obs is given by (|3.4p and (|3.5p for the Gold and 
SNLS datasets respectively. 

The steps we followed for the minimization of (j3.6p for 
the Gol d and SN LS datasets are described in detail in 
Refs [9i[ini[ioi- The form oiH{z) used in (1^, (1^ 
is obtained by using the CPL parameterization (|2.18p and 
leads to 

H^{z) = Hl[no;u{l + zf + 

+ (l-l]o,n)(l + ^)'''+"°+"^'ewf] (3.9) 

The best fit form of w{z) for the Gold and first year SNLS 
datasets is shown in Fig . 2 for a prior of ilom — 0.24 along 
with the Icr errors jllOj l (shaded region). Even though 
the two fits are consistent at the la level it is clear that 
the Gold sample mildly favors a crossing of the PDL at 
z ~ 0.25 while no such trend appears for the SNLS data. 
This trend difference could be attributed to three factors 

• Statistical errors 

• Systematic errors 

• Different redshift ranges (Gold < z < 1.7, SNLS 
< z < 1). 

The third possibility has been excluded in Ref. |lOl| 
where a truncated version of the Gold sample with (0 < 
z < 1) was found to show even stronger trend for crossing 
the PDL than the full Gold sample. Given this appar- 
ent difference in trends between the Gold and the SNLS 
datasets the following question arises: 'Which (if any) 
of the two trends do other cosmological probes of dark 
energy favor?' 

B. The CMB shift parameter 

A particularly accurate and deep geometrical probe of 
dark energy is the angular scale of the sound horizon at 
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the last scattering surface as encoded in the location if'^ 
of the first peak of the CMB temperature perturbation 
spectrum. By measuring the angular scale Of^ ~ V^f"^ 
of the last scattering sound horizon and calculating its co- 
moving scale rs{zdec) independently from CMB physics, 
the sound horizon angular diameter distance dA{zrec) can 
be obtained assuming flatness as 



dA^Zrec) 



1 



1 



dz 
Wz) 



(3.10) 



thus providing a useful constraint on H{z). To eliminate 
the model dependence involved in the calculation of the 
sound horizon scale an (approximately) model indepen- 
dent parameter can be defined by dividing the measured 
angle 6^^ ~ 1/ by the corresponding angle 0^'^ of a 
reference model. This parameter is known as the 'shift 
parameter'. 

The shift parameter is defined as [lOi, [lo^, [lOl [ml 



R = 



'TT 



TT 



(3.11) 



where if^ is the temperature perturbation CMB spec- 
trum multipole of the first acoustic peak. In the defini- 
tion of R, if"^ corresponds to the model (with fixed il-om, 
^06 and h) characterized by the shift parameter and l-^'"'^ 
a reference fiat SCDM model (fio„j = 1) with the same 
uj,n = flomh^, ujb = flobh^ as the original model. 

The location if^ of the first acoustic peak can be con- 
nected with the angular diameter distance dA to the last 
scattering surface and with the sound hori zon r^ at the 
last scattering surface (z = z^ec) as follows |lllLlll2| |: 



where 



ll = lAidA,rs)il - ^(uJrajUJb)) 



lA{dA,rs) = -p™ 



dA^Zrec) 



(3.12) 

(3.13) 
(3.14) 



and the phase shift parameter tp{uJm,^b) — 0.27 
depends weakly on cosmological parameters. The sound 
horizon r^ and the angular diameter distance depend 
on the Hubble expansion history (at early and late times 
respectively) as follows: 

Cs{a)da 



am{a) 
Cs{a)da 

n 



1/2 



-1/2 



(3.15) 



where Cs{a) is the sound velocity which at decoupling 
(a = flrec) is 



Sp 

Tp 



1 



1 



31 



(3.16) 



3000 
2000 

; 1500 

1000 




FIG. 3: Derivation oftfie sfiift parameter: The 3-year WMAP 
binned data along with a theoretical flat ACDM model (curn ~ 
0.14, = 0.022, /i = 0.72) and a SCDM model with {uj„ = 
0.14, tJt = 0.022, /i = 0.37). For the ACDM the first peak is 
at ij'^ = 220 while for the SCDM model at = 246. 



and 



dA{z) 



zsm 



dz' 



For the reference model SCDM we have 
lo 



(3.17) 



Ho 



K„(i + z)3 + f7;(i + z)4]i/2 



(f7; + 1)1/2 -(a;,, + f^;)i/^" 



(3.18) 



Using equations (|3.12p - (|3.18p it is easy to show that 



I'TT 



R 



4(1-^') 



ir uii - ^) 

r^d^4(4ec) _ 2 q{ni,arec) 

r's dA{Zrec) f^J/f Xi^) 



where 



Xiz) 



Hodz' 



(3.19) 



(3.20) 



and we have assumed flatness and used the fact that 
i) = -0' (since ujh = uj'b and LOm = w^) and arec{oJb,u!m) = 
^'reci^'b^^'m) [13 ■ The main advantage of using the 
shift parameter instead of just the value of ij"^ is that 
it is weakly dependent on parameters other than ge- 
ometry {H(z)). The dependence on other parame- 
ters enters through ilgi^ and through q{Q'^,arec) = 

q{uJr/h''^,arec{'-Om,tUb}) whcrC UJr = Q,rh? ~ 2.47- 10~^, 



and 



1 



1048(1 
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0.00124W, 



-0.738 



){l + giuj9^) where gi = 0.0783w 



-0.238 



(1- 



39.5wj^-^63)-i and 52 = 0.560(1 + 21. Notice 
that the weak dependence of the shift parameter on Wm, 
cjf, through q(ujm,i^b) has not been exphcitly demon- 
strated in previous studies [sH]- Instead, the shift pa- 
rameter is usually expressed as 



R' = 



^Om Xi^rec) 



— 1/2 

or equivalently as i? = ^Q^^xi^rec) thus omitting the 
correction factor q{u!m,^b)- For fixed r^om — 0.27, fib = 
0.043 and 0.5 < h < 0.9, q{ujrn,^b) remains in the range 
0.964 < q < 0.968 and therefore the dependence on the 
Hubble parameter introduced by q is very weak. 

The main disadvantage of the shift parameter is that 
in order to evaluate it using its definition (j3.1ip for a 
given CMB perturbation spectrum we need not only the 
location of the first peak if^ but also the location of 
the reference model (SCDM) peak. For an accurate re- 
sult, the latter requires running a Boltzmann code like 



(3.21) 



CAMB hlj with n'r 



Om 



1, ujL 



0J„ 



ojb and 



evaluating the location of the first peak 1-^^'^ . This pro- 
cedure is illustrated in Fig. 3 where we show the 3-year 
WMAP data along with a theoretical fit obtained for 
{Lu,n,LOb,h) = (0.14,0.022,0.72) assuming flat ACDM. 
The first peak for the best fit model is obtained at 
If^ = 220 ± 0.8 while for the reference fiat SCDM 
model with (w™,u;b,/i) = (0.14,0.022,0.37) (leading to 
= 1) we have l'-^^ = 246 (see Fig. 3). We thus 
obtain 



R 



'TT 



iTT 



1.123 ±0.03 



(3.22) 



Using now ((HT^ and ([X^ with q = 0.965 we find 

i?= rij/fx(^rec) = 1.71 ±0.05 (3.23) 

This result is consistent but slightly different from that 
of Ref. [35| which used MCMC chains to obtain R from 
WMAP3 as 



R= 1.70 ±0.03 



(3.24) 



In what follows we adopt the result of eq. p.24p to obtain 
the best fit w{z). 

Even though the shift parameter incorporates most 
of the geometrical information of the CMB spectrum 
the following question can be raised: 'How much of the 
CMB spectrum geometrical information is left out of the 
shift parameter?' To answer this question we proceed in 
two steps: First we consider the CPL parameterization 
and construct CMB spectra (using a modified version of 
CAMB[113|]) varying H{z) with fixed R and other CMB 
parameters (see Fig. 4 and Table |T]). The constructed 
CMB spectra are practically identical for all values of I 
except low / {I ~ 50) where there are minor differences 




FIG. 4: Theoretical CMB spectra for the CPL parameteri- 
zation for five models that have fixed R — 1.70 and other 
parameters (see text and Table |T]|. Models (1,2,3) with iden- 
tical Qom and R but different ^1(2:) have practically identical 
spectrum for all / (identical ISW effect). The ISW effect for 
fixed R changes only due to variations of Qom (models 4, 5)! 



due to the ISW effect [401, [llj. These CMB spectra dif- 
ferences can not be observationally distinguished due to 
the large cosmic variance errors that dominate at low I. 

The temperature anisotropy due to the ISW effect is 
given by an integral over the time variation of the gravi- 
tational potential traversed by the photons coming from 
the last scattering surface f40| 

ei(fc,77o) = (2/ + 1) / dr,e-^ 2$(fc,7;) ji{k{i^a - n)) 

(3.25) 

where the overdot indicates conformal time derivative, 
Q = ^ and the integral is taken over the line to the last 
scattering surface. Due to the integral over the oscillat- 
ing Bessel functions j; , the dominant contribution to the 
CMB anisotropy is on large scales (low I). The gravita- 
tional potential <I>(fc, 77) (metric fluctuation of FRW back- 
ground) may be shown (using the Einstein equations) to 
obey a Poisson equation which in a flat universe is of the 
form 



where 5p is the density perturbation 
dominated universe v-cTm §£. 



(3.26) 
In a flat, matter 



we 



have ^ 

p 



I =^ Sp ^ a 

which implies (due to p.26p ) that $ = and there is 
no ISW effect. However, in the presence of dark energy, 
the evolution of perturbations is delayed at late times 
< . ^^^"^ therefore $ 7^ leading to a non- 
zero ISW effect. The ISW effect therefore measures the 
growth rate of perturbations at late times or its deviation 

from the flat matter dominated rate ^ ^ a. 

p 

An interesting feature of Fig. 4 is that if we keep R 
and ilom flxed while varying H{z), the ISW effect appears 
to be unchanged! In contrast, if we keep R fixed while 
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energy models for two reasons: 

1. The shift parameter has a very low relative error 
(few percent) 

2. The integral over H{z) extends to high redshift 
{zrec — 1089). Therefore, minor modifications of 
H{z) can be very 'costly' in exactly due to this 
term. 

Another similar term is obtained from the imprint of 
the last scattering somid horizon on the large scale struc- 
ture correlation function. 



FIG. 5: The growth of perturbations D(a) for the five models 
of Fig. 4 along with the corresponding forms of w{a). Those 
models (1,2,3) which have identical ISW effect (same Qom and 
R but different H(z) — ™(z)) also have very similar growth of 
perturbations even though the dark energy density evolution 
is quite different. 



varying flom and H{z) {px{z)), the ISW effect changes 
indicating that perhaps there is not much more infor- 
mation in the ISW effect about the dark energy density 
evolution than that included in the shift parameter R. 
This is also demonstrated in Fig. 5 where we show the 
growth of perturbations D{a) for the five models of Fig. 
4. Clearly, those models which have identical ISW effect 
(same ilom and R but different H{z)) also have very sim- 
ilar growth of perturbations even though the dark energy 
density evolution is quite different! 

It should be stressed that whatever information about 
dark energy is encoded in the ISW effect, this information 
can not be extracted directly from the temperature per- 
turbation CMB spectrum (due to the large cosmic vari- 
ance errors at low I) but only using alternative methods 
like CMB cross corre latio n |40l| with large scale structure 
or CMB polarizationflTsj. 

The measurement of the shift parameter (equation 
(I3.23p ) allows us to add an important term to the 
of equation p.6p obtained with Snia. This term is of the 
form 

2 .o ^ (j?(»oin,wo,u;i)-1.70)^ 

XcMB{^Ora,WQ,Wi) = -^-^ (3.27) 

This term is very important in constraining evolving dark 



TABLE I: Parameters for the 5 models of Fig. 4. 



Models 




TOO 




h 


Model 1 


0.27 


-0.8 


-0.123 


0.725 


Model 2 


0.27 


-0.9 


0.196 


0.725 


Model 3 


0.27 


-1.0 


0.497 


0.725 


Model 4 


0.14 


-1.7 


-1.061 


0.998 


Model 5 


0.50 


-0.3 


-0.064 


0.533 



C. The Baryon Acoustic Peak 

The large scale correlation function measured from the 
luminous red galaxies spectroscopic sample of the SDSS 
(Sloan Digital Sky Survey) [1] includes a clear peak at 
about lOQh-^Mpc. This peak was identified with the 
expanding spherical wave of baryonic perturbations orig- 
inating from acoustic oscillations at recombination. The 
comoving scale of this shell at recombination is about 
ISOMpc in radius. The identification of the comov- 
ing scale where the correlation peak was observed with 
the acoustic horizon at recombination requires the ap- 
propriate form of the cosmological model H{z) in con- 
verting from the observed angles-redshift to correlation 
function distances. An accurate determination of the 
best fit H{z) would therefore correspond to using a gen- 
eral model or parameterization of H[z) to convert from 
the observed angles-redshifts to correlation function dis- 
tances and then vary these parameters until the observed 
Baryon Acoustic Oscillation (BAO) peak coincides with 
the expected value from CMB physics. 

A simplified version of this approach was used in Refs. 
(3, [lil]. It involves the use of a fiducial ACDM H{z) 
model to construct the correlation function from the ob- 
served angles-redshifts and identify the distance scale 
fpeak where the BAO peak appears. Comparing this 
value of Tpeak with the expected value r'^^^^, from CMB 
physics would require a scale shift by a factor a: r'^^^f. = 
ctrpeak- This factor a can be approximated as the ratio 
of the required characteristic distance scale Dv{zbao) 
of the survey with mean redshift zbao (obtained with 
the correct H{z)) over the distance scale Dy^^^'^ {zbao) 
corresponding to the fiducial ACDM model 

Dv{zbao) 



peak 
^peak 



nKCDM 



{zbao) 



(3.28) 



The characteristic distance scale of the redshift sur- 
vey with typical redshift Zs can be connected to H{z) 
as follows: Consider a spherical shell (Fig. 6) of comov- 
ing radius R. Let z be the redshift corresponding to the 
sphere (points C and D) as viewed by an observer at O 
and Az the redshift difference between A and B. 

If we measure the angular scale A0, the redshifts z and 
Az, then given a cosmological model H{z), the comoving 
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FIG. 6: A spherical shell in redshift space. 



scales CD and AB may be evaluated from the flat FRW 
metric as [ll7j 



CD = xM = M 



dt _ Ae r dz 

lu a. Ho Jq E{z) 
where E{z) = H{z)/Ho and 

M 1 Az 

AB = Ax = 



(3.29) 



a Ho E{z) 



(3.30) 



Generalizing the above to a deformed sphere [AB ^ CD) 
we may associate an approximate single scale to the 
sphere as [3] 



Dv{z) 



( CD 



-I 1/3 



1 



{AB) 



dz 



Ho Jo E{z)J HaE{z) 



1 A2 



1/3 



(3.31) 



where CD/ Ad is the comoving angular diameter dis- 
tance. We also want this scale to be representative of 
the whole sample (fiducial scale for SDSS zbao — 0.35) 
so we allow the sphere to extend from the observer to 
[A ^ O) to z ~ 0.35 {zbao — 0.35) and therefore 
Az ~ Zbao — 0.35. For low redshift samples (like 
Zbao = 0.35), the dilation scale Dv{z) can be assumed 
to encode all the information required for converting red- 
shift to distance. By demanding that the scale of the 
peak of the correlation function coincides with the acous- 
tic scale (which is predicted by simple physics (see ec 
(|3.15p ) and depends on ^omh"^) it has been shown 
that 

ZJy (0.35) = 1370 ± UMpc (3.32) 
for ^omh^ = 0.130 ± 0.010. A dimensionless and inde- 



pendent of Ho version of the dilation scale Dy is 



A = Dv{0.i^) 



\/^OmHf 

Zbao 



= nli^E{zBAo)-'/'[ 

= 0.469 ±0.017 



l2/3 _ 



ZBAO Jo E{z') 



(3.33) 



This parameter has been used e xtens ively in constraining 
dark energy models (see eg [t^ Ill8l |) and it will also be 
used in what follows. 

The measurement of the parameter A (equation (|3.33p ) 
permits the addition of one more term in the of equa- 
tions p.6p . p.27p to be minimized with respect to H{z) 
model parameters. This term is 



XBAoi^Om,Wo,Wi) 



(A(Oom,wo,wi)- 0.469)2 



0.0172 



(3.34) 



D. The Cluster Baryon Gas Mass Fraction 

The baryon gas mass fraction for a range of redshifts 
can also be used to constrain cosmological models H{z). 
The basic assumption underlying this met hod is that the 
baryon gas mass fraction in clusters @, Ill9l | : 



fgas 



Mb- 



gas 



(3.35) 



is constant, independent of redshift and is related to the 
global fraction of the universe in a simple way. This 
relation may be written as 



3 

^On 



Mb 



(l + a) 



Mb 



Mt, 



(l + a) fgas (3.36) 



where 6 is a bias factor suggesting that the baryon frac- 
tion in clusters is slightly lower than for the universe as a 
whole. Also 1 -I- a is a factor taking into account the fact 
that the total baryonic mass in clusters consists of both 
X-ray gas and optically luminous baryonic mass (stars), 
the latter being proportional to the former with propor- 
tionality constant a ~ 0.19\/7i [6]. Assuming now that 
the hot gas in a cluster follows a spherically symmetric 
isothermal model i.e. 



"-eo(l + -^) 



_^-3/3/2 
2 - 



(3.37) 



where /? is a constant, is the electron number density 
and rieo, ?'c are the central electron density and the core 
radius respectively, it is straightforward to show that [36l . 
[T20t 

Mgasi< R) = B{T,,Rlr,,p)rl/^Lx{< R)^'^ (3.38) 



In eq. (|3.38p Tg is the electron temperature, Lx{< R) 
is the bolometric luminosity within radius R and i? is a 
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TABLE II: Cluster data from Ref. i6i 





J g as,i 




0.078 


0.189 


0.011 


0.088 


0.184 


0.011 


0.143 


0.167 


0.019 


0.188 


0.169 


0.011 


0.206 


0.180 


0.015 


0.208 


0.137 


0.018 


0.240 


0.163 


0.009 


0.252 


0.164 


0.012 


0.288 


0.149 


0.017 


0.313 


0.169 


0.010 


0.314 


0.175 


0.023 


0.324 


0.177 


0.018 


0.345 


0.173 


0.019 


0.352 


0.189 


0.025 


0.363 


0.159 


0.017 


0.391 


0.159 


0.024 


0.399 


0.177 


0.017 


0.450 


0.155 


0.019 


0.451 


0.137 


0.009 


0.461 


0.129 


0.019 


0.461 


0.156 


0.034 


0.494 


0.094 


0.025 


0.539 


0.135 


0.011 


0.686 


0.155 


0.018 


0.782 


0.100 


0.016 


0.892 


0.114 


0.021 



constant independent of cosmological parameters. The 
quantities Tc, Lx{< R) are obtained from the observed 
angular core radius 9c and the observed apparent lumi- 
nosity within 6 {lx{< d)) as 



Lx{< R) 
R 



4TTdLizflxi< 

OcdA{z) 
OdAiz) 



(3.39) 
(3.40) 
(3.41) 



where dL^z) and d^iz) are the luminosity and diameter 
distances which depend on the cosmological model H{z) 
as 



dL{z) = {l-VzfdA{z) = {l + z) 



dz' 
W) 



(3.42) 



in the context of a flat cosmology. From eqs (|3.38p - (|3.42p 
we obtain 



Mgas{< R) = c{ec,e,ix.z)dA{zf'^ 



(3.43) 



where all cosmological model dependence is encoded in 
d^(z) while C depends on observables characterizing each 
cluster (6'c, lx{< 0),Te,z). 

Similarly, the total cluster mass within R Mtot{< R) 
may be obtained assuming that the intracluster mass is 
in hydrostatic equilibrium as 



Mtot{< R)^DR = D9 dA{z) 



(3.44) 



where D is independent of cosmological model parame- 
ters. From eqs p. 431) and (|3.44p we find that 

where Q depends only on individual cluster observables. 
Assume now that the quantities Qi {i = 1,...,N) have 
been obtained for an observed sample of N clusters at 
redshifts Zi. Using eqs. (|3.36p and (|3.45p we have 



= (1 + a)fgas{z^) = (1 + a)Q,dA (2.)'^' (3-46) 



(3.47) 



Define now 

/JL'^^*'^(z,)^Q.4^^*^(2.)^/^ 

where d^^^^''{zi) is the angular diameter distance cor- 
responding to SCDM (flat r^om = 1) used as a refer- 
ence model. Solving (|3.47p for Qi and substituting in eq. 
((Oel) we find 



jSCDM 
"A 



(2.) 



dA{Zz) 



3/2 



(3.48) 



rSCDMi N ^ 
J gas 1^2 j — 

Since now fg^s^'^^ is known observationally (see eq. 



1 + a f^On 

cSCDM 
gas 

p.47|) ) we can use eq. (|3.48l) to find the best fit dA{zi) 
and therefore the corresponding best fit cosmological 
model. We use the 26 cluster data for fgas^^^{zi) pub- 
lished in Ref. 6] and minimize Xcbf (Cluster Baryon 
Fraction) defined as 

26 / fSCDM/^ \ ^ f ^2 
2 \ \J qas V-^jy Jgas,i) 
XcBF[^''am,Wo,Wi) = 2 

i=l fgas.i 

(3.49) 

where fgas^^'^izi) is given by eq. (|3.48p and the ob- 
served fgas,i, f^gas i uscd in our analysis are shown in 
Table [III We treat as a nuisance parameter and 

we marginalize over it as follows: Define A = tt— tt^ 
then 



SCDM 
gas 



\ J / \ / — "J qas i^i) 
\ dA{Zi) J 

(3.50) 

and eq. (|3.49p after expanding with respect to A becomes 



xIbf = A'^ -2XB + C 



(3.51) 



where 



N 



A 



B 



C 



E 

i=l 
N 

E 

i=l 
N 

i—l J a 



fSCDMf^ \2 
J gas 



J qas K'^i) Jgas 



(3.52) 
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Equation (I3.5ip has a minimum for A = i?/^ at 

xIbf = C- By a (3.53) 

which is independent of A. The Xcbf P-53P is the 

one actually used in our analysis. It should be pointed 
out however that the relative errors of these datapoints 
are more than 10% and therefore variations of H{z) with 
given r^om prior affect the value of x%bf much less than 
Snia, CMB and BAO data. 



Our goal is to connect the galaxy overdensity in redshift 
space 5s = with the corresponding galaxy over- 

density dg = ^''"^^^ in real space. 

The number of galaxies in a spatial region is the same, 
independent of the coordinate system i.e. 



ns{xs)d^Xs = n{x)(f'x 



(3.60) 



or 



E. Linear Growth Rate at z=0.15 



n(l + (53) =n(l + (5g)^^ 



(3.61) 



The 2dF galaxy redshift survey (2dFGRS) has mea- 
sured the two point correlation function at an effective 
redshift of Zs = 0.15. This correlation function is af- 
fected by systematic differences between redshift space 
and real space measurements due to the peculiar veloc- 
ities of galaxies. Such distortions are expressed through 
the redshift distortion parameter (3 which connects the 
power spectrum in redshift space Ps{k) with the true 
galaxy power spectrum Pg{k) as 



(3.54) 



where ^ = cos 9 and 9 is the angle between k and the line 
of sight. The parameter (3 can be observationally deter- 
mined from the observed power spectrum (or its Fourier 
transform, the correlation function) in redshift space and 
may be shown to be connected to the growth factor D{a) 
at the effective redshift of the power spectrum as [121] 



where 



^4 



dD(a)/da 
g = a- 



D{a) 



(3.55) 



(3.56) 



where n is the average number density of galaxies. As- 
suming low velocities and focusing on modes fca; > > 1 we 
obtain from eqs (|3.59p and (|3.6ip 



(3.62) 



By Fourier transforming eq. (|3.62p and expressing the 
peculiar velocity in terms of the mass overdensity 5 using 
the continuity equation 



a— — h ikv — 
dt 



which implies [1211 



v{k, a) 



igaH6{k, a)k 



we find (for low redshifts z << 1) 

5s(k) = 5g{k) + g^iH{k) 



(3.63) 



(3.64) 



(3.65) 



where ^ = cos 9 = k ■ x and g is defined by eq (|3.56p . We 
now use the bias factor b = to write eq. (|3.65p as 



and b is the bias factor {b = which connects the over- 
density in galaxies Sg with the matter overdensity S. Th e 
derivation of eq. p.55p may be sketched as follows 
The radial position of a galaxy with low redshift z and 
no peculiar velocity is approximated by 



X{z) = jr (3-57) 

-no 

and therefore its redshift space position vector is 

2; 

Xs = — —{sin 9 COS (j), sin 9 COS if), COS 9) (3.58) 
Ho 

Due to gravitational effects however, galaxies have pe- 
culiar velocities v which affect the redshift z as z = 
Hqx+v-x where x is the line of sight direction. Therefore, 
the true comoving distance x to the galaxy is connected 
to the redshift inferred distance as 



z 

i^o 



V ■ X 



(3.59) 



l{k) = {l + ^fi^)Sg{k) (3.66) 

which leads to eqs p.54p and (I3.55p . 

The parameter P — j may be measured from redshift 
surveys by measuring Ps{k) and expanding both sides of 
eq. p.54p in Legendre polynomials 7'/(^)[12l]: 

Ps{k) = Pi°)(fc)7'o(M) + PP{k)V2{f^) + Pi'Hk)V4{f^) = 
= (1 + Pii^fPg{k) - [(1 + + i/32)7'o(M) + 

+ {\P +\p^)'P2{^i) + ^P^V^{^i)]Pg{k) (3.67) 

Therefore one way to measure (3 is to use the measured 
quadrupole Pi^^(fc) and monopole Pi°^(fc) as 



1R2 



pr{k) i + i/3 + i/3^ 



(3.68) 
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which leads to the value of p and (if the bias is deter- 
mined) to the value of g at the effective redshift of the 
redshift survey. 

The correlation function (the Fourier transform of the 
galaxy power spectrum) can also be used instead of the 
power spectrum to obtain [3 using a very similar method 
as the one described above. Such a method was used 
along with others in Ref. [1] to find 

/3 = 0.49 ±0.09 (3.69) 

at the effective redshift of z = 0.15 of the 2dF redshift 
survey. This result may now be combined with the linear 
bias parameter 

6=1.04 ±0.11 (3.70) 

obtained from the skewness induced in the bispectrum 
of the 2dFRGS[37t by linear biasing to find the growth 
factor g at z = 0.15 as 

= a^^U=o.i5 ^^'•/S = 0.51 ±0.11 (3.71) 

D{a) 

Such a result can in principle be used in combination 
with eq. (|1.14p with f{k,a) = 1 to constrain H{z) in 
the context of general relativity or, if H{z) is determined 
by other observations, to determine the function f{k, a) 
in eq. p.l4p thus providing a robust test of alternative 
gravity theories. 

Under the assumption of general relativity, the mea- 
surement of the Perturbations Growth Rate (PGR) g{z = 
0.15) can be used to add one more term to as 

2 /o ^ (5r(f7om, Wo, M^l) - 0.51)2 

where g{Qoai,wo,'Wi) is obtained by solving equation 
(|1.14p with f{k, a) — 1 and initial conditions D{a) ~ a 
for a ~ 0. 

Notice however that the large relative error of this 
datapoint (about 25%) makes the contribution of Xlss 
rather insensitive to variations of H{z). If however future 
observations reduce the relative error of this measure- 
ment it is potentially very important as it is a dynamical 
test capable of distinguishing between general relativity 
and extended gravity theories along the lines discussed 
in section I. 



F. Cosmological Data and PDL Crossing 

In order to identify the trends encoded in the current 
cosmological data with respect to dark energy evolution 
and in particular PDL crossing we group the data in three 
categories: 

1. SnIA Gold sample 

2. SnIA SNLS 



3. Other Dark Energy Probes (ODEP) that include 
the CMB, BAG, Gluster Baryon Fraction (CBF) 
and Perturbations Growth Rate (PGR) at z = 0.15 

The corresponding to be minimized for each category 
are 

• XGoid(^Om,u'o, wi) (equation ([^ ) 

• x|7VLs(^Om,t«o,Wi) (equation ((3^ ) 

• Using equations dS^T]), ([glM]) . ^^M\f and ((X7^ we 
have 

XoDEP = XCMB + XbAO + XCBF + XpGR (3.73) 

In order to identify the dependence of the resulting 
best fits on the fiom prior used we have not marginalized 
over ilom- Instead we have fixed Qom and considered two 
cases (rJom — 0.2 and ilom = 0.3). The range between 
the two cases includes the current best fit value of f2om 
based on WMAP and SDSS which is[12^] Qom = 0.24 ± 
0.02. The errors of the best fit wjz) were obtained using 
the covariance matrix (see eg |l23l Il24l |). The best fit 
form of w{z) for each dataset category is shown in Fig. 
7 for both fiom = 0.2 and fiom = 0.3. The corresponding 
contours in the wq — wi parameter space is shown in 
Fig. 8. 

The following comments can be made based on our 
results summarized in Figs. 7 and 8: 

• The Gold dataset mildly favors dynamically evolv- 
ing dark energy crossing the PDL at z ~ 0.2 over 
the cosmological constant while the SNLS does not. 

• Dark energy probes other than Snia mildly favor 

crossing of the PDL for low values of fiom {^Om ^ 
0.25) while for Oom — 0.3 this trend is significantly 
reduced. 

• The best probes of dark energy with currently ex- 
isting data, other than SnIa are the GMB shift pa- 
rameter and the BAG peak. 

The last point emerges by plotting the third column of 
Figs. 7, 8 utihzing only GMB and BAG data. Such plots 
are practically unchanged compared to those of Figs. 7, 
8. 

An important issue that needs to be addressed before 
closing this section is the issue of systematics introduced 
by using a particular parametrization to fit the cosmolog- 
ical expansion history. There is no doubt that the best 
fit forms to a given dataset may differ significantly be- 
tween different parametrizations. However, the particu- 
lar property of PDL crossing at best fit has been shown to 
be a robust feature for several H{z) — w{z) parametriza- 
tions in the context of the Gold dataset provided that 
these parametrizations allow for PDL crossing. This was 
demonstrated in Ref [9^ (Fig. 1) where we compared 
the best fit forms of a wide range parametrizations to 
the Gold dataset. 
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0.6 1.2 1.8 0.6 1.2 1.8 0.6 1.2 1.8 




0.6 1.2 1.8 0.6 1.2 1.8 0.6 1.2 1.8 



FIG. 7: The best fit form of w{z) for each dataset category for both Qom ~ 0.2 and ilom = 0.3 along with the la errors (shaded 
region). The categories are: Gold dataset (column 1), SNLS (column 2) and Other Dark Energy Probes (ODEP column 3). 




Wo Wo Wo 



FIG. 8: The 68% and 95% confidence contours in the wo — wi parameter space for each dataset category for both Qom = 0.2 
and f2om = 0.3. Notice that for the SNLS dataset ACDM is within the la region. 
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In that paper it was demonstrated that even though 
the high redshift properties at best fit were different for 
each parametrization, the low redshift properties and 
in particular the PDL crossing which appears at z ~ 
0.2 are very robust and consistent among the different 
parametrizations. This robustness can also be checked 
in the context of Other Dark Energy Probes. We have 
fitted three different H{z) parametrizations to the ODEP 
along the lines of Ref. |99|. We found that the best fit 
forms for each parametrization indicate crossing of the 
PDL for low Ooin but in this case there is a spread in the 
redshift range of crossing in the range 0.4 < z < 1.8. In 
view of the very broad redshift covering (0 < 2 < 1089) 
of the dominant probe of this class which is the CMB 
shift parameter, the above spread in the crossing range 
is not a significant indicator of systematics. 

Another source of systematics that deserves attention 
is the robustness of our results shown in Figs 7 and 8 
with respect to the theory of gravity considered. Mod- 
ified gravity theories are expected to affect dynamical 
tests of dark energy (which depend on the way the met- 
ric couples to energy momentum tensor) but not geo- 
metrical tests which detect directly the evolution of the 
FRW-metric. Tests based on standard candles (Snia) 
and standard rulers (CMB shift parameter, BAO peak 
and X-ray gas mass fraction) are geometric tests and we 
do not expect modified gravity theories to have a direct 
effect on them. For example even though the CMB spec- 
trum will vary in modified gravity theories (in particular 
its ISW part), its geometric features like the first peak lo- 
cation depend only on the background metric and act as 
a standard ruler measuring directly the integral of II(z). 

On the other hand the test based on the growth rate 
of perturbations at z = 0.15 would give a different form 
for H{z) in the context of modified gravity theories (see 
Fig. 1). This effect however would hardly modify our 
results of Figs. 7, 8 because the la error of the pertur- 
bations growth contributions to is much larger that 
the geometrical tests contributions. In fact we have found 
that even if we completely ignore the dynamical test con- 
tribution to Figs. 7, 8 these figures remain practically 
unchanged. We should stress however that there is an 
indirect effect of modified gravity theories on SnIa stan- 
dard candles. The possibility of a time varying Newtons 
constant would induce a variation of the Chandrasekhar 
mass and lead to a variation of the S nIa a bsolute luminos- 
ity. This effect was studied in Ref. (l25j for an evolving 
Newton's constant consistent with nucleosynthesis and 
solar system bounds. It was found that the effect is neg- 
ligible for current SnIa data but may have to be taken 
into account for future data with SNAP level accuracy. 

It is therefore clear that even though currently avail- 
able dark energy probes agree on the fact that the dark 
energy equation of state w(z) is close to —1 there is no 
universal trend with respect to the evolution properties 
of dark energy and the existing trends appear to depend 
on the value of r^om prior considered. 



IV. CONCLUSION-OUTLOOK 

We have pointed out the importance of identifying 
observational signatures that could distinguish between 
dark energy and extended gravity theories as possible ori- 
gins of the accelerating expansion of the universe. The 
crossing of the PDL w — —1 can be identified as such 
a signature based on a geometrical dark energy probe. 
We have noted that the simplest theoretically motivated 
class of theories that is consistent with such crossing is 
the class of extended gravity theories. A representative 
example of this class of theories includes scalar-tensor 
theories. 

Even though there are currently no clear observational 
indications for a PDL crossing, there are probes of ac- 
celerating expansion that mildly favor such crossing over 
a constant equation of state. These probes include the 
Gold SnIa dataset for any reasonable prior for fJom as 
well as the combination of other than SnIa probes (CMB, 
BAO, Clusters Baryon Fraction and growth rate of per- 
turbations) with a low Jlom prior (f)om ^ 0.25). On the 
other hand a non-evolving w ~ — 1 appears to be favored 
by the SNLS first year data. 

Another (perhaps more robust) signature which could 
rule out dark energy in favor of extended gravity theories 
is obtained by utilizing the measurement of the growth 
rate of perturbations at various redshifts. This test uti- 
lizes a 'dynamical' (as opposed to 'geometrical') dark en- 
ergy probe. We have demonstrated how can this mea- 
surement be used to distinguish between the two classes 
of theories. We have also pointed out that current 2dF- 
GRS data are consistent with ACDM in the context of 
general relativity and require no extended gravity theory. 

We have also discussed the ISW effect as another 
potentially interesting dynamical probe of dark energy 
which is sensitive to the growth rate of perturbations 
at recent redshifts. We have demonstrated that in the 
context of general relativity and for fixed shift parame- 
ter R and fixed ilom the ISW effect on the CMB power 
spectrum and the growth factor of perturbations D{a) 
appear to be insensitive to modifications of the dark en- 
ergy equation of state. Degeneracies of the CMB power 
spectrum with fixed shift parameter R have been dis- 
cussed previously |114| but in those cases fiom was varied 
simultaneously with a w{z) ~ wq — const to keep R con- 
stant. Thus, the large scale part of the CMB spectrum 
was seen to vary in those studies due to the variation of 
D{a) (ISW effect). The degeneracy observed here can 
be used as an additional observational discriminator be- 
tween general relativity and extended gravity theories. 
For example if the ISW effect is found to differ from the 
expected form obtained by the measured shift parameter 
R and fJom then this could be viewed as an indication of 
extensions of general relativity. A more quantitative form 
of this argument is an interesting potential extension of 
this work. Modifications of the ISW effect are expected 
for example in theories predicting a time varying New- 
ton's constant. In particular in Ref. [l26| | it was shown 
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that in the DGP model at late times, perturbations enter 
a DGP regime in which the effective value of Newton's 
constant increases as the background density diminishes. 
This leads to a suppression of the ISW effect, bringing 
DGP gravity into slightly better agreement with WMAP 
data than conventional LCDM. 

The PDL crossing viewed as a geometrical signature 
of extended gravity assumes the non-existence of phan- 
tom degrees of freedom which is a reasonable assump- 
tion in view of the theoretical problems of such degrees 
of freedom. On the other hand the dynamical signature 
(growth rate of perturbations) assumes that the dark en- 
ergy perturbations can be ignored on sub-horizon scales 
and can not mimic the effects of modified gravity theories 
expressed through the function f{k, a) (equation (jl.l4|) ). 
These assumptions make it important to identify further 
types of observational signatures that can be used in com- 
bination with the above, providing more robust tests that 
could distinguish between the two classes of theories. 

Numerical Analysis: Our numerical analysis 
was performed using a Mathematica code implement- 
ing fitting of multiple dark energy probes and a 
modified version of CAMB allowing for a dynamical 
w{z) using the CPL parametrization. All the codes 
and data along with detailed instructions are available at 



http: / /leandros. physics. uoi.gr/pdl-cross/pdl-cross. htm 
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